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Abstract 

Causal Dynamical Triangulations is a background independent approach 
to quantum gravity. In this paper we introduce a phenomenological trans- 
fer matrix model, where at each time step a reduced set of quantum states 
is used. The states are solely characterized by the discretized spatial vol- 
ume. Using Monte Carlo simulations we determine the effective transfer 
matrix elements and extract the effective action for the scale factor. In 
this framework no degrees of freedom are frozen, however, the obtained 
action agrees with the minisuperspace model. 

1 Introduction 

The model of Causal Dynamical Triangulations (CDT) was proposed some years 
ago by J. Ambj0rn, J. Jurkiewicz and R. Loll with the aim of defining a lattice 
formulation of quantum gravity from first principles [IJ [2] . The foundation of 
this model is the formalism of path-integrals applied to quantize a theory of 
gravitation. The quantum gravity path integral is regularized by discretizing 
the spacetime geometry g with piecewise linear manifold T. The building blocks 
of four dimensional CDT are four-simplices, which properly glued along their 
faces form a simplicial manifold. 

An important assumption of CDT is the causality condition. As a conse- 
quence of the original Lorentzian signature of spacetime, only causal geometries 
should contribute to the integral. We will consider globally hyperbolic pscudo- 
Riemannian manifolds which allow introducing a global proper-time foliation. 
The leaves of the foliation are spatial three-dimensional Cauchy surfaces called 
slices. Because topology changes of the spatial slices are often associated with 
causality violation, we forbid the topology of the leaves to alter in time. For 
simplicity, we chose the spatial slices to have a fixed topology of a three-sphere, 
and establish periodic boundary conditions in the time direction. Therefor, the 
spacetime topology is M = S 1 x S 3 . Spatial slices of a triangulation are enu- 
merated by a discrete time coordinate t. To each vertex of the triangulation 
such time coordinate is assigned, bringing on a distinction between space-like 
links of length a s and time- like links of length a t . Because each simplex contains 
vertices lying in two consecutive spatial slices, there are two kinds of simpliccs: 
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first of a type {4, 1} with four vertices lying in one spatial slice and one in the 
neighboring slice, and second of a type {3, 2} with three vertices lying in one 
spatial slice and two in the adjacent slice. The Wick rotation is performed by 
the analytic continuation to imaginary lengths of the time- like links a t — > ia t . 
The regularized partition function Z is now written as a sum over causal trian- 
gulations T, 

z= f v[ g y sEH ^ -+ J2 e ~ s[T] - (!) 
j r 

The Einstein-Hilbert action S EH [g] = J dt fd s Xy/=g(R - 2 A) evaluated 
on a simplicial manifold T composed of N4 simplices, among them N41 being 
of type {4, 1}, and N vertices, gives the discrete Regge action, 

S[T] = -K N + K 4 N 4 + A (N 41 - 6JV ), (2) 

where Kq, K 4 and A are bare coupling constants, and naively they are functions 
of G, A and at,a s - 

We applied Monte Carlo techniques, and using the Regge action @, mea- 
sured expectation values of observables within the CDT framework. The sim- 
plest observable is the scale factor a(t), or more conveniently the three- volume 
nt defined as the number of tetrahedra building slice t. 

For a certain range of the coupling constants, a typical configuration is bell- 
shaped, with the average volume profile (n t ) oc cos 3 (</-B). The emerged back- 
ground geometry behaves like a well defined four-dimensional manifold and is 
perfectly consistent with a Euclidean dc Sitter universe, the classical vacuum 
solution of a spatially homogeneous and isotropic minisuperspace model [4]. In 
earlier work we have shown [3] that the discretized minisuperspace action, 

S[n t ] = ~J2 ( (nt+1 T ? + ^ - X A < 
r t \ n t+i + n t ) 

describes well not only the measured {nt) but also the fluctuations 

C tt ' = ((n t - (n t ))(n t , - (n t ,))). (4) 

The effective action (j3]) couples only adjacent slices. Such form suggests that 
there exists an effective transfer matrix labeled only by the scale factor. 

2 The transfer matrix 

The model of Causal Dynamical Triangulations is completely determined by a 
transfer matrix Ai labeled by three-dimensional triangulations r. The matrix 
element (ti|.M|t2) denotes the transition amplitude in one time step between 
states corresponding to triangulations n and r-x- It is given by the sum over 
all four-dimensional triangulations T of a slab, with boundary triangulations t\ 
and T2, 

(ti\M\t 2 )= e ~ S[n - 
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The transfer matrix M. depends both on the entropy factor, which counts num- 
ber of triangulations T connecting the boundaries in one time step, and the 
Regge action S[T]- The partition function ([T]) corresponding to T time steps is 
then expressed in terms of the matrix M., 

Z = Y / ^ Sm =^M T . (5) 
r 

The probability of finding a configuration with T spatial slices given by three- 
dimensional triangulations n , t 2 , . . . , tt is 

p( T )(n, . . . ,r T ) = |(ri|M|T 2 )(75|M|r 3 ) . • . (T T \M\n). (6) 

We used partition function ((SJ) in Monte Carlo simulations. The measure- 
ments performed so far, have been concentrated on the measurement of the 
three-volume n t . The probability p( T >(n\, . . . ,tit) of finding a configuration 
with spatial volumes n\, n 2 , • • • , n-r is given by a proper sum of partial prob- 
abilities ((BJ). Let T 3 (n) denote the subset of three-dimensional triangulations 
which are build of exactly n three-simplices. We use the projection operator 
p{n) = \n)(n\ on the subspace spanned by T$(n), 

p{n) = \n){n\= £ |r>(r|. (7) 

reT 3 (n) 

to express the probability (rii, . . . , Ut), 

pV\n x ,...,n T ) = ^T r [\n 1 )(n 1 \M\n 2 )(n 2 \M\n 3 )...(n T \M}. (8) 

In ([5]) it is misleading to think of the aggregated "state" \n) as a normalized 
sum of the vectors |t), r 6 23(71). Such a vector would again be a single vector 
located in the space spanned by the |r)'s. It is more appropriate to interpret 
the "state" associated with n as arising from a classical uniform probability 
distribution of states |r) and in this way to treat p(n) as the associated density 
operator. 

As mentioned in the Introduction, the form of the effective action ([3]) ob- 
tained from the covariance matrix §4§ suggests that there exists an effective 
transfer matrix (n\M \m) whose elements are labeled by the three- volumes and 
that it is possible to effectively decompose observed distributions {ri\, . . . , rir) 
into a product 

P^ T \n u ...,n T ) = i(ni|M|na>(na|M|n3)-"(nr|M|ni>. (9) 

The effective transfer matrix M depends only on the coupling constants Kq, 
A and K4 but not on the number of slices T. In analogy to ([6]), the elements 
of the effective transfer matrix correspond to transition amplitudes in one time 
step between states of a given three-volume. 
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Figure 1: Left: The empirical transfer matrix for range 1200 < n t < 1600. 
Right : The difference between the empirical and theoretical matrices disappears 
in the numerical noise. 



3 Measurements 

In the following, we will assume that we can work with an effective transfer 
matrix (n\M\m) and will show that equation ([9]) provides a very good approxi- 
mation of measured data [5]. 

For simplicity, let us define the two-point function, 

P iT \n t ,n t+At ) = ^(n t \M At \n t+At )(n t+At \M T ~ At n t ). (10) 

by summing (jH]) over all three-volumes except for times t and t + At. The 
simplest way to measure the matrix elements (n\M\m), up to a normalization, 
is to consider T = 2, 



(n\M\m) cx y^p( 2 )(m = n,?i2 = m). 



The effective transfer matrix elements can be measured in various ways. In 
particular, as follows from ([TUf for T = 3,4, we have 

{n\M\m) cx 1 ' (11) 

y "W(ni = n, 71.3 = in) 

We tested, that the elements (n\M\m) measured in different ways completely 
agreed up to numerical noise, supporting validity of equation (|9]). For technical 
reasons, most measurements were performed using expression (|lll) . 

The coupling constant A'4 in ([2]) plays a role of a cosmological constant. To 
correctly perform simulations, we have to approach with K4 very close to its 
critical value K^ rlt . To efficiently probe desired range of the three- volume, we 
added to the Regge action @ a quadratic term to fix n t around n vo i, 

S -> S + e^J(n t - n vo i) 2 . 
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Figure 2: Left: (n|M|c — n) plotted as a function of n for various c (dots). 
Gaussian fits are drawn with a line. Right: the coefficient k(c) in the kinetic 
term as a function of c = n + m, (different colors denote different ranges) and 
a linear fit k{n + m) = T ■ (n + rn) (red line). 



Because it is consistent with the decomposition ©, its effect can be easily 
canceled. For technical reasons, we measured the transfer matrix M separately 
for few overlapping ranges of the three- volume. 



4 The effective action 

The effective action obtained from the covariance matrix (U) is directly related 
to the effective transfer matrix M. The minisuperspace action §3§ suggests that 
the effective transfer matrix given by 



(n\M\m) =We L v J (12) 

is a good approximation in the bulk where n t is large 0. Further, we will mea- 
sure the empirical transfer matrix elements (n\M\m), extract the parameters T, 
fj, and A, and check that (|T2l is indeed a good approximation of the data. The 
measured effective transfer matrix M, for range 1200 < n t < 1600, is presented 
in Fig. [1] (left graph) . The right graph shows the difference between the mea- 
sured matrix M and the best fit (fT2]) . Indeed, the difference disappears in the 
numerical noise proving that the approximation (|12[) is very good. 

The measurements presented in this paper were performed for coupling con- 
stants K = 2.2, A = 0.6 and K A = 0.922. 



4.1 The kinetic term 

To get a better estimation of the parameters associated with the effective action 
^ and (fT2")l , we first try to fit only to the parameters of the kinetic term which 

x We slightly modified the form of the potential term. Such parametrization is more con- 
venient to extract the parameters of the action. 
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Figure 3: log(n|A/|n) of the scaled transfer matrix (dots, different colors denote 
different ranges) compared with the fit of the potential term — £ e // (red line, 
which stops at n = 400). 



is the dominating term from a numerical point of view. We do that by keeping 
the sum of the entries, i.e. n + to, fixed such that the potential term is not 
changing. In this way we determine T with high accuracy. The matrix elements 
for constant n + m = c show the expected Gaussian dependence on n (see left 
Fig. ®, 



(n\M\m) = (n\M\c - n) = Af(c) exp 



(2n - c) 



(13) 



where the terms in the effective action which only depend on c are included in 
the normalization. 

We expect the denominator of the kinetic term k(c) to behave like k(n+m) = 
r ■ (n + to). As shown in the right graph in Fig. [5] this is indeed true and the 
parameter T is constant in the whole range of the three- volumes. The best linear 
fit gives r = 26.07 ± 0.05. This result is consistent with the values obtained 
from the fits for separate ranges of n t . 

4.2 The potential term 

The potential part of the effective Lagrangian L e f / may be extracted from the 
diagonal elements of the transfer matrix, 

L e ff(n,n) = — log(n|M|ra) + const = i (jin 1 ^ - An) . (14) 

For technical reasons, we measured the transfer matrix M separately for few 
different ranges of the three- volume. Because, the normalization is not uniquely 
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defined, in order to merge the effective Lagrangian, the constant in ([T4"]) has 
to be properly adjusted. The measured merged effective Lagrangian is shown 
in Fig. [3] The colors denote different ranges for which the transfer matrix was 
measured. Fig. [3] presents also the fit of form (fT4")l . In the bulk region, where nt 
is large enough, the theoretical expectation (THJ fits very well. The measured 
values arc fi = 16.5±0.2 and A = 0.049±0.001, where we took T = 26.07. Again, 
this result is consistent with the values obtained from the fits for separate ranges 
of n t . 



5 Conclusions 

The model of Causal Dynamical Triangulations comes with a transfer matrix 
(ti|A^|t2). The measured distributions of the three- volumes n t , e.g. P( T \n t , n t +At), 
have an exact definition in terms of the full transfer matrix M. and the density 
matrix |n)(n|. The actual data coming from Monte Carlo simulations seem to 
allow for a much simpler description in terms of an effective transfer matrix M , 
labeled by abstract vectors \n) referring only to the three-volume. The effective 
transfer matrix M allows to directly measure the effective action /Sfnt]. An 
important advantage of the present method, since number of slices T is small, 
is much faster measurement of the transfer matrix compared to the covariance 
matrix, which was used previously to extract the effective action. Basically over 
the whole range of nt the effective transfer matrix elements can be represented 
as 



(n\M\m) =Ne r 

with high accuracy. This result is fully consistent with the reduced minisuper- 
space action (j3|), although in CDT we do not freeze any degrees of freedom. 

An issue not addressed in this article, is the problem of small three- volumes. 
For small n t we do not observe a Gaussian distribution of the three-volume nt 
around the mean value (n t ). Because of strong discretization effects, the prob- 
ability distributions, and consequently the effective transfer matrix elements, 
split into three families [SJ. Despite different nature, after the smoothing pro- 
cedure, the effective action for small volumes is basically the same as for large 
volumes, with a small modification in the potential [5J. It might be interpreted 
as possible curvature corrections, however, we are not able to measure it accu- 
rately in a discretization independent way. 
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